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FLUID LIMIT OF A HEAVILY LOADED EDF QUEUE WITH 
IMPATIENT CUSTOMERS 

L. DECREUSEFOND AND P. MOYAL 


Abstract. In this paper, we present the fluid limit of an heavily loaded Earli¬ 
est Deadline First queue with impatient customers, represented by a measure¬ 
valued process keeping track of residual time-credits of lost and waiting cus¬ 
tomers. This fluid limit is the solution of an integrated transport equation. 
We then use this fluid limit to derive fluid approximations of the processes 
counting the number of waiting and already lost customers. 


1. Introduction 

Queueing theory is a keystone of the development of current telecommunications 
systems. Engineers now aim to guarantee the grade of service customers are enti¬ 
tled to receive according to their contracts with the carrier. One way to meet this 
objective is to schedule requests according to their “importance”. Of the utmost in¬ 
terest, are the audio and video traffic flows, which are subject to severe transmission 
delay constraints. In particular, some requests can be thought as “impatient” since 
it may be better to discard some packets which would eventually arrive too late in 
order to favor some other packets which still can meet their delay requirements. 

Another very active branch of queueing theory is, nowadays, devoted to the 
analysis of call-centers where customers are impatient: they tolerate to wait up to 
a certain limit upon which they depart from the queueing line hence are considered 
as lost both for the queueing system and for the called-service provider. It is then 
crucial to develop service disciplines which ensure a maximal number of served 
customers by controlling waiting times keeping them within impatience bounds. 
These disciplines are commonly referred to as real-time service disciplines. 

In real-time queuing theory, each customer is not only identified by his arrival 
time and service duration but also by a deadline. This means that a customer has 
a given period of time (his time credit, i.e., the remaining time before his deadline) 
during which he should enter the service booth. This time credit decreases at unit 
rate as time goes on. If it expires before the customer enters service, the customer 
is either lost and the discipline is said to be hard, or he is kept in the waiting 
line and the discipline is said to be soft. The discipline we address here is the so- 
called Earliest-Deadline-First discipline in its “hard” version: the customer having 
the smallest time credit is served first and whenever the credit-time of a customer 
expires before it is served, this customer is lost. 

To find which service discipline is best, one can compare them within a static 
scenario, i.e., customers to be served are all present at initial time and no new 
customer enter the system, service duration and impatience of each customer are 
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all known at the beginning; or in dynamic environments, i.e., customers arrive 
randomly, their service duration and impatience are only known stochastically. In 
both settings, it appears that the so-called Earliest-Deadline-First (EDF for short) 
discipline is optimal. It is known for a while |Der M that EDF discipline is optimal 
for the static approach: if any (real-time) service discipline can serve the customers 
of a given scenario without loss then EDF also does. Within random environments, 
it has been proven in |PT88| and generalized in |Moy05| that EDF discipline ensures 
the least possible failure probability, i.e., the least number customers lost by missing 
their deadline. 

Yet, apart from the notable exception of deterministic deadlines for which EDF 
discipline reduces to the FIFO service policy with impatience, no closed form of the 
loss probability is known. The only satisfying quantitative approach so far consists 
in numerically assessing the loss probability for an EDF system with Markov-chain 
approximations tHXD88l[PK9lllND^|PT88l . 

When no simple tractable object can describe a queueing system, one wants to 
identify its “mean behavior”. One hopes that a Markovian process characterizing 
the system, when suitably normalized, can be approximated by a fluid limit that 
is, a deterministic continuous function of the time. Then the fluid limit describes 
the general behavior of the considered process. Numerous queueing systems have 
already been investigated this way (see, for instance, |RohfinirRor()7] for a pure delay 
system, [DLSfll] for a soft EDF queue, |OPWf)l] for a queue run under a processor 
sharing service discipline). For instance, Lt being the amount of customers at time 
t in an M/M/1 queue with parameters A and fi, one proves that the sequence of 
processes defined by = 1 and for all t > 0, := n~^Lnt tends 

in distribution to L := ((1 -|- (A — /i)t)^, t € K''") and that {-\/n 
converges in distribution to a diffusion process. The fluid approximation of the 
system presents the same first order characteristics as the “real” system: it fills in 
at velocity A and empties at velocity the congestion reaches 0 to the condition 
A < /i (this is Loynes’s stability condition) after a time A — /i (mean duration of a 
busy period). 

We want to obtain the same type of information for an M/M/1 queueing system 
with impatient customers. In this case, it is easily seen that the process 
which counts the number of customers in the system is no longer Markovian. Indeed, 
the value of Xt+h not only depends on Xt, but also on all the time credits of all the 
Xt customers present in the queue at time t. Therefore we describe the system by 
the point measure-valued process (t't)(>o whose unit of mass are the time credits 
of all the customers waiting in the queue, or already discarded. 

Formally, it is rather straightforward in our case to write down the infinitesimal 
generator of the Markov process {yt)t>o^ see Theorem |2l below. It is made of four 
terms, all but one are standard and represent the evolution of the process when an 
arrival, a departure or nothing occurs during an infinitesimal time period. The nat¬ 
ural but unusual term is the term due to the continuous decreasing of the residual 
deadlines at unit rate as time goes on. This term involves a “spatial derivative” of 
the measure v, a notion which can only be rigorously defined within the framework 
of distributions. Because of this term, the fluid limit equation (see is the inte¬ 
grated version of a partial differential equation rather than an ordinary differential 
equation as it is the rule in the previously studied queueing systems. Thus, the 
famous Gronwall’s Lemma is of no use here. Fortunately, the partial differential 
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equation which pops up, known as transport equation, is simple enough to have 
a closed form solution - see Theorem ^ Thanks to that, we can then proceed as 
usual to show the strong convergence of the renormalized process to the fluid limit. 

This paper is organized as follows. After some preliminaries, we deflne and solve 
the integrated transport equation in the space of tempered distributions. In Section 
0 we establish that the above described process (ut)j>Q, is a weak Feller Markov 
process and give its inflnitesimal generator. In Sectional we prove the fluid limit 
theorem. The last section is devoted to applications to the EDF driven queue with 
deterministic initial time credits, and to a pure delay system. 


2. Preliminaries 

We denote by V},, respectively Cq and Cb the set of real-valued functions deflned 
on R which are bounded, right-continuous with left-limits (rcll for short), respec¬ 
tively continuous vanishing at inflnity and bounded continuous. The space Vb is 
equipped with the Skorokhod topology and Co and C& with the topology of the uni¬ 
form convergence. The space of bounded differentiable functions from R to itself is 
denoted by C^ and for (j) £ Cl, |j (j) ||oo:= sup(|^(a;)| -I- |( 5 i'(a;)|). For all f £ T>b and 

all a; S R, we denote by Txf, the function Txf(-) '■= /(■ — x). 

The Schwartz space, denoted by S, is the space of infinitely differentiable func¬ 
tions, equipped with the topology deflned by the semi-norms: 

I (j) |a,h:= sup I I, a G N, 5 e N. 

a,gK dx° 

Its topological dual, the space of tempered distributions, is denoted by S', and the 
duality product is classically denoted (/i, (j)). The Fourier transform on S is deflned 
by := (27r)“^/^ Jg^e~''^^(f>{x) dx and the Fourier transform is deflned on S' by 
the duality relation {Jl,(j)) = {fi, (j)). 

The set of finite positive measures on R is denoted by and A4p is the set of 
finite counting measures on R. The space Af^ is embedded with the weak topol¬ 
ogy, a(A4'j ,Cb), for which is Polish (we write (/i, /) = j f dfx for /i G Alj and 
/ G I?h). We also denote for all a: G R and all v £ AiJ, TxV the measure satisfying 
for all Borel set B, Txv{B) ■= v {B — x) . Let Co(AlJr,R), be the set of continuous 
functions from to R, vanishing at inflnity, endpwed with the topology of the 
sup norm. Let 0 < T < oo, for E a Polish space, we denote C ([0,T],i5), respec¬ 
tively V the Polish space (for its usual strong topology) of continuous, 

respectively rcll, functions from [0,T] to E. 


(E) 


3. The integrated transport equation 
The transport equation on C^(R x R+,R) with unknown u{x,t) is deflned as: 
dtu = —bdxU -I- / in R X (0, oo), 
rt = h at R X {f = 0}, 

where 6 is a real number, / is a function of (R x R-|-,R), and h £ (R, R). It 

is well known (see |Eva,98j l that (E) admits a unique solution given for all x,t by: 

(1) u{x,t) = h{x — tb) + f f{x + {s — t)b,s)ds. 

Jo 





4 


L. DECREUSEFOND AND P. MOYAL 


Let US define the following extension of the transport equation: 

Definition 1. Let T > 0, K G S', igt)t>o ^ ([0, T],5') such that go = 0 and 

b be a real number. The process {gt)t>Q satisfies the integrated transport equation 
E(K,g,h) on V ([0, T], 5') if for all (f & S', and for all t G [0, T]: 

E(if, g, b) {r)t, 4>) = {K,4>) -b f {gsA') ds + {gt, 4>). 

Theorem 1. The integrated transport equation \E(K, g, bjjj admits a unique solu¬ 
tion (Lt)t>o ^ ([0 iE],5'), satisfying for all (j) G S and for all t G [0,T]; 


(2) 


{Lt, (j)) = {K, Tbt (j)) + {gt A) -b / {gs, n^t-s) <^') ds. 

Jo 


Proof. Let L and M be two solutio ns of (E(iL, g, b) I and let iV = L — M. For all 
t € [0, r], it follows from E(jL, g, b) that for all </) G 5: 

l(N,A) = -b{NtA'). 

Denoting for all C G •= this can be rewritten: 

d 


dt 


{NtA) = -b{Nt,iffi) = -b{NtA<f) = -bANtA). 


Solving the latter differential equation yields for all (/> G 5 and alH G [0,r]: 

(NtA) = {Noe^^^cfi) = {No,e^^A) = (No,^^) = 0 , 

hence for all t G [0,T], Nt = 0. Therefore, there is at most one solution to 
< E(jL,g,b)t . 

The process (Lt)(>Q defined by Q belongs to V {[0,T],S') and for (j) G S, we 
have: 

b f {LsA') ds 
Jo 

= b {K,Tbs4>') dsb {gsA')ds-b^ j / {gr,Tb(s-rA') dr ds. 

Jo Jo Jo Jo 

Since dt{C,TbA) = -b{C,,nA')^ we get: 

bj {LsA')ds = -J ^{{K,Tbs(t))) ds + b j {gsA')ds 

+ bj J ^{{gr,Tb(^s-rA')) dsdr 

= -{K,TbA) + (KA) +b f {gs,Tb(^t-sA') ds = -{LtA) + {KA) + {gtA), 

Jo 


The process (Lt)i>o thus satisfies (E(iL, g, 6) I. 


□ 
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4. The profile process 

Following Barrer’s notation |Ba,r57j . we throughout this paper consider a queue¬ 
ing system with impatient customers M/M/l/l+GI-EDF: 

• customers arrive at times The process defined for all t by 

■= ^ ^{Ti<t} 

ieN* 

is a Poisson process of intensity A > 0, 

• a first sequence of marks {cilieN*; sequence of service durations re¬ 
quested by the customers, is i.i.d. with the distribution of tr which an 
exponentially distributed with parameter /i > 0 random variable, 

• the customers are impatient: i.e., the i-th customer leaves the system, and 
is lost forever, when he doesn’t reach the service booth before his specific 
deadline, Ti + Di. In other words, he is initially labelled with a random 
variable referred to as his patience, or initial time credit, Di. The marks 

independent and identically distributed with the distribution 
of D, an almost-surely non-negative and integrable random variable. The 
time credits of the customers decrease continuously with time, at velocity 
one (in time units). Provided that the i-th customer entered the system 
before t {Ti < t), but did not reach the service booth before t, we denote 
Di{t) the residual time credit at t of this customer, i.e., the residual time 
before his possible elimination. Therefore: 

Di{t) = Di - {t- Ti), 

and Di{t) < 0 means that the *-th customer has been lost, reaching his 
patience before t before entering the service, 

• there is 1 non idling server and a buffer of infinite capacity, 

• the service discipline is EDF (i.e.. Earliest Deadline First): when completing 
a service, the server deals with the customer whose residual time credit is 
the smallest among all the customers in the buffer, if any. This service then 
proceeds until completion, without any interruption. 

Let us finally define the following performance processes: 

Xt := Number of customers in the system (buffer + service booth) at t , 

Qt := {Xt — 5')+ = Number of customers in the buffer at t , 

St := Number of customers served up to time t , 

Pt := Number of customers lost up to time t , 

At time t, provided that the buffer is non-empty, denote for z = 1,...., Qt, Ri{t) the 
z-th residual time credit of a customer in the buffer at t, ranked in the increasing 
order: 

Ri{t) < i?2(t) < .... < RqAt)- 

Provided that at least one customer has been lost at t {Pt ^ 0), for z = 1,. ,Pt, 

denote R-i{t), the z-th residual time credit among the customers lost up to t in the 
decreasing order: 


R-pA^) < R-Pt-iA) < < R-i{t)- 
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The time credit profile of the system at t is the following measure: 

Qt Pt 

i—1 i=l 

where Sx is the Dirac mass at x. Provided that the buffer is non-empty at t, we 
denote for all i = 1, Qt, 

U{vt) ■■= Ri{t), 

the i-th point of ut (in the increasing order) on the positive half-line. 

The service discipline can be represented as follows : when the server completes 
a service (say at time s), he first deals with the customer whose time credit is given 
at this time by: 

ti{Rs) = Rl{s) > 0, 

provided that Qs ^ 0. The customer corresponding to the atom fi(Ps), being chosen 
by the server, leaves the buffer: the corresponding atom is erased from the 

point measure vt for all f > s (this customer won’t ever reappear in the buffer, since 
the service discipline is non-preemptive). 

By profile process of the queue, we mean, the process {i't)t>o credit 

profiles at t. This process is fully characterized by its initial value uq, the real 
numbers A > 0 and ^ > 0 and the non negative integrable random variable D. 
This process will consequently be referred to as the profile process associated to 
{pq, X, fi,D). The dynamics of the profile process can be depicted as follows. The 
atoms are translated continuously towards left at velocity 1, at the arrival time R, 
an atom is added to the measure vt^ at Di the initial time credit of the arriving 
customer, and at an end of service Ti , an atom disappear from the measure Uf, at 
Figure ^ shows a typical path of the profile process. Note, that the buffer 
congestion and loss processes can be deduced from the profile process by writing 
for all f > 0: 

Qt = (l^ijllR^), Pt = (j^tjlR-), 

since the waiting, resp. already lost, customers at t are those who have positive, 
resp. non positive, time credits at t. 


5. Markov property 


Denote for all t, At the remaining time before the next arrival after t, and for 
all t such that Xt > 0, Ft the remaining time before the next end of service after t. 
For all t,h> 0: 


vt+h = < 


ThVt 

Ph^t ^ti{ut)—h 
PhVt + 


if At > h and Ft > h, or Qt = 0, 
if At > h, Ft < h and Qt > 0, 
if At < h and Ft > h ot Qt = 0 
and the customer arrives at tk 
affected with the initial time credit dk, 


the more complex events (several arrivals, several ends of service, or arrivals and 
ends of service) between t and t + h being of probability o{h). This dynamics 

shows in particular that {vt)t>o ^ ^ (l®’ ’ ®i^ce ((z/t, 4’),t> 0) belongs to 

D([0,oo),M) for all (j) G C/,. We finally define the filtration: 


Ft := CT (i^s(R), s <t,B G iB(R)). 
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Beginning of service between t and t + h 


0 Ri { t ) R 2 (.^) 

^^H-H--h- 


Ri{t+h) 


0 ' 

H- 




Lost customer 


Ri{t+2h) R3{t+2h) 

f R2(t+2h) 


Arrival between t and t + 2h 


t + h 


t + 2h 


Figure 1. Dynamics of the profile process. 


Theorem 2. The profile proeess associated to (uq, X,pi,D) is a weak Feller 

process with respect to (.?^t)t>o j whose infinitesimal generator is given by: 


(3) AF{iy) = ^im - (a + 

+ nF{iy- ,5t,(^)) l{^(R^)>o} + A J F{R + 5d) dPoid), 

for all F in the domain of A: 

F{A) :=Co(Ad/,M) n “IliRi —- exists. 

Proof. For all t,h>0 and all bounded measurable function F : A4^ ^ R: 


( 4 ) 


E 


[F{Rt+h)\At] =(l - (A + /rl|^^(„^)>o|) h^F{ThRt) 

+ fihF [rhVt — ThSti(ut)) l{yt(R* )>o} 

+ \h J F {ThRt + ThSd) dPoid) + o{h) 

=■■ ThF{Rt). 


Thus, according to |Da,w93j . p.l8, [vt)t>o ^ weak homogeneous Markov process, 
whose transition function is given by (T/j, h > 0). For F € Cq , it is easily 

seen from 1^, that ThF G Cq for all h > 0. Since A4'^ embedded with 
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the weak topology is locally compact separable, it routinely follows that (Pt)t>o is 
a weak Feller process whose infinitesimal generator of v is given by □ 

Corollary 1. For all (j) G Cl, the process defined for all t > 0 by: 

(5) M^(t) = {vt, (f) - (pq, (!>)-[ {vsA') ds 

Jo 

is an rcll Ft-martingale such that M,j,(t) G for all t > 0. Its increasing process 
is given for all t > 0 by: 

(6) < >t= p f (1)^ {ti (vs)) 1{^,(r;)>o} ds + AfE [(/i^(iA)] . 

'I 0 

Proof. Let (p G Cl. Define the mapping II^ : K for all v by: 

n^iv) := {v,4>). 

Since 

( 7 ) -{v,(t>'), 

we have for all u G Mpi 

AU,f,{L') = -{v, (j)') - pP (fi(p)) l{^(R^)>o} + AE [0(D)] . 

Furhtermore, 

^(n0(T/,p)^ -n0(p)^)= -2(p, 0)(p, 0'), 
and hence for all v G Mpi 

Anl{v) = 2(p, 0)yfn0(p) + pffi {h{v)) 1{^(R. )>o} + AE [0^(D)] . 

From Dynkin’s lemma |EK8fil|D^ n H, it follows that and the process defined 
for all t by 

- Al{vo) - [ Anl{v,) ds 
Jo 

are J^t-local martingales. This entails that <M^ >t=< {v.,4>) >t for all t > 0, and 
thus, Ito’s integration by parts formula yields to: 

(pt,0)^ = (po,0)^ + 2 / {vs,(t>) dM^{t) 

Jo 

+ 2 (Ps,0)Ain,^(Ps)ds+<(p.,0)>t . 

Jo 

Hence, for all f > 0: 

2 / {vs,(t>) dM^{t)+<{v,,(t))>t 
Jo 

= N^t) +p f (p^ 1{^,(r;)>o} ds + AtE [02(D)] , 

'J 0 

and by identifying the finite variation processes, we obtain 0 . □ 
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6. Fluid limit 


For all n G N*, denote by Uq , a measure on and define (u")j>g, the profile 
process of the M/M/l/l+GI-EDF queue whose initial state is represented by the 
profile Vq, whose arrival process (-/V")(>o is Poisson of intensity A" >, where the 
customers request service durations are exponentially distributed of mean expecta¬ 
tion and have initial time credits i.i.d. with the distribution of iA". The 

process (u")(>g is in other words the profile process associated to (ug , A", lA"). 

Also denote > 0), the associated filtration, 

rg":= inf{t>0,ur(M;) =0}, 

the first time when the buffer is empty, and 

Wo := inf {< > 0,ti(z/t") = 0}, 

the first time of loss of the system. We also define as previously the performance 
processes of the n-th system: (Ai")j>g, (QDoo! given for all thy Q'^ = {v'^, 1r. ), 
given by P/" = (u^", 1 r_). 

According to Theorem for all (f) the process defined for alH > 0 by: 


(8) M;{t) = (ur, 4>) - (Ug", 4>) - /Vr, ^') ds 

Jo 

+ ds - A"fE [</>(P")] 

is an rcll iF"-martingale, such that M^{t) G for all t, and whose increasing 
process is given for alH > 0 by: 

(9) < M; >,= (ti(ur)) l{,n(R;)> 0 } ds + A"tE [</>2(P")] . 


We normalize the process (u")j>g in time, space and weight the following way: for 
all Borel set B and for all t, define 


^t{B) 


KtinB) 

n 


where 


nB := {nx, x G P} . 


The first positive atom of u" is therefore given by: 




p {<t) 

n 


We also denote > 0) := > 0), the associated filtration, 

fg" := inf {t > 0, ur(R+) = O} = ^Tg", 


uiq := inf {t > 0, ti(u") = 0} = —ujq 

n 

and normalize the arrival process as well as the performance processes of the n-th 
system the corresponding way, i.e., for all t > 0, 




/\rn _ vn 

-\rn ^nt /~\n 

—A t — A Qt 


Q n pn 

nt pn _ -^nt 

t ' ‘ t ' 




10 


L. DECREUSEFOND AND P. MOYAL 


For all t > 0, Q" Pp can thus be recovered by: 

(10) = 


(11) pr = (^r,iK-)- 

Let (j) G Cl and i/’"(.) = /^- As easily seen from (|H1) and 0, the process 

defined for all t by 

(12) M-{t) := M;„(nf) = {i?-, 4>) - ct>) - [\p:, <t>')ds 


+ m" / .^(ti(pr))l|,n(R.)>o}ds-A"tE 


pn 

n 


is a t/”-martingale of I? ([0, oo), R), such that G Lp' for all t. Its increasing 

process is given for all t by: 


/■* A” 

< Mil (^i(^r')) 1 { p "( r :;)> o } ds H-fE 




We now define the set of hypothesis under which we will prove a law of large 
numbers for the sequence of processes p”. 


Hypothesis 1. 


• There exists two real numbers ^ > 0 and X > ^ such that: 


(13) 


A"^A, 

n —>00 


(14) 


d —^ P- 

n —>oo 


• For all e > 0, there exists > 0 such that for all n G N*, 

(15) P [(Po”, 1) > nM,] < e. 

• There exists a measure Pq of j\4ji such that for all / G Vb: 

V 


{(^>0./)} 


nGN* 


'4J)- 


• There exist an integrable and almost surely non-negative r.v. D such that: 

V 


D, 


E 


D^- 


E [D] . 


Proposition 1. Assume that Hypothesis^ holds. Then, \ (r'”)(>Q t is tight in 

t - J kGN* 

I?([0,oo),M+). 

Proof. According to Jakubowski’s criterion |Daw93j . it suffices to show that: 

(1) For all (f) G Cl, the sequence {((^"1 </’))t>ol tight in I? ([0, 00 ), R), 

I — JnGN* 

( 2 ) For all T > 0 and 0 < 77 < 1, there exists a compact subset Kt^t/ of A4ji 
such that 


liminfP [p" G KprfVt G [0,T]] >1-77. 











FLUID LIMIT OF A HEAVILY LOADED EDF QUEUE 


11 


In order to prove the first condition, let us fix € Cl and T > 0. Remarking that 
for all 5 > 0 , 

(i^:,l)< V" + (Po",l), 

Equation 111 ‘211 yields for all u < u < T: 

(16) < [ I</>')I rfs + ds 

J U j u 


+ A”E 


on 

n 


<\v-u\ II <l)' Iloo + k - u\ II Iloo {Vq, 1) 

+ \v-u\\\ (j) Iloo (Im" - mI + |A” - A|) + |u - u 


(A + fi) 


Let e > 0 and 77 > 0. First, let 

From Markov’s inequality: 

P 


:= 


er] 


30 II (j)' Iloo XT' 


sup |u - u| II cj)' Iloo > 7 

,V<.T ^ \v — u\<Si ^ 


< 


rj 

AT5(5i II (/)'Iloo , |A-A"|T55 i II (/)' 


g |A-A"|T5<5i II cj)' 
6 


V 


and thus with there exists A^i > 0 such that for all n> Ni, 


sup |u - m| II (j)' Iloo > 


il,v<.T , |v —u|<(5i 


< £ £ ^ £ 

“ 6 6 3' 


Let now: 


<52 := 


According to 111 511 . for all n G N* 


5M,/3 II </.' 


sup |u-u| II (j}' Iloo (1^0 > 1 ) > r 

U,V<.T , \v — u\‘<52 ^ 


< P 


(^oM)> 


V 


= p [(PoM) > m,/3] <-. 


° " 5<52 II (j)' 

According to assumptions ltT3|l and itTHl , there exists N 2 , such that for all n> N 2 , 
for all u,v <T, 

|u-u| ||</,||^(|7,"-7,| + |A"-A|)<|, 

and letting 


<5s := 


V 


6 


(A + 77 ) ’ 
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sup III - m| II Iloo (A + Ai) < (^3 II (/' Iloo (A + At) < -. 

u,d<T , |d —u|<(53 ^ 

Now, let ^ > 0. Apply successively Markov’s and Doob’s inequalities: 


sup|M^(t)| > ^ 


t<T 


< [< MJ >t] 




,,n fT 

I a 2 


— ^ (^(^”))1{pj‘(r;)>o}^« + 


<1^/^ 

-e^ 


—TE 

n 

n yn 

n n 


D" 

n 


T —>0. 

n— 9-00 


For all n G N*, > 0^ being a rcll process on [0, T], one can apply the stan¬ 
dard convergence criterion |Robnnj . from which it follows that | ^ o) | 


n^N* 


converges in distribution to the null process. This sequence is in particular tight in 
V ([0, T], K.): there exists ^4 > 0 and N 3 > 0 such that for all n > N^: 


sup \M;iv)-M;{u)\>{ 

<i,v<T, \v—u\<64 ^ 


< 

- 3 


Finally, in view of the previous inequalities and CEl, there exists 5 > 0 and TV G N 
such that for all n > A^: 


sup - {iy^,(j))\>r] 

u,v<T , |v—n |<(5 


(17) P 

On the other hand, let 

Ue ■■= Me II (j) 

Assumption m implies that for all n € N*, 


< £. 


(18) 


P[|(^o",<^)|>a,]<P 


</.||oo Me II ^ 


< e. 


With itT^ and itTTIl we can apply the standard tightness criterion of real valued 
processes (see for instance [Rohnf)| 'l: for all T > 0, < y 4 ')) t>Q r tight in 

V ([0, T], R): it is tight in V ([0,00), R). 

We now prove the second tightness condition (compact containment). Let us 
first apply |GPW01j. Lemma A.2.: under hypothesis^ we have the following weak 
law of large numbers: 

JV" 


1 




(AiEp(S)]) inUdO.T], 


new 

for any (j) G Cb- In particular, this yields for any 0 < I <T: 

. ^n(t + i) 


sup - ^ <^(A") > 2A1E [0(.D)] 

telo,T-i]n 


(19) P 

Taking I = T and ^ = 1 in the last expression yields: 


0. 
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Denote /(.), the identity on R. Taking I = T and (j) = I in also leads to: 


N" 


Let 


- E A" > 2ArE [D] 

Mt = max|2Ar+ (P^,l),2ArE [D] + (up, J)| + 1. 


sup max{(u”, 1 k+), (u", J1 r+)}> Mt 
te[o,r] 


We have: 

(20) P 

< P [(u", 1) + > 2AT + {D*,l) + 1] 

1 " _ 

P (9^,1) + - V > 2ATE [D] + {9*J) + 1 


2=1 


< P [(U-, 1) > (p*, 1) + 1] + P [N^ > 2AT] 


+ P[(Po",/)>(P*,/) + l]+P 


> 2ATE [D] 


i=l 


0 . 


Let us now define, for all T > 0, 0 < 17 < 1, the set 

I^T,rj ■= |c G -Adj ; max{(^, 1r+), (C,/1r+)}< Mt , {(, l(_oo,-r]) = o|- 

Since (C,/1 r+) < Mr, this implies that for all y > 0, ( ([y,oo)) < Mr/y, and thus 

lim sup C([y, oo ))=0 , lim sup C ((- 00 , y])) = 0 , 
y^^CelCT.r, y^-^C&ICT.r, 

which implies that ICtti C Af is relatively compact IkM . Now, since up to time 
T no lost customer can have a residual time credit less than —T, 

sup(P”,l(_oo-T]) = 0. 

t<T 

This, together with piji implies that: 


lim inf P 


u" G ICt, 7 j, for all t G [0,T] 


>1 — 77 . 


Kt, 7 ) being the closure of /Ct,?;, we found a compact subset Kt,,; C Af ^ such that: 


lim inf P 

n—^oo 


G KT,r„ for all t G [0,T] 


>1 — 77 . 


□ 


Theorem 3 (Fluid limit theorem for M/M/l/l+GI-EDF queues). Assume that 
Hypothesis^ holds, that there exists T > 0 such that: 


(21) 


p K <T]-^o, 
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and that there exists a deterministic element (Pt)t>o o/C ([ 0 , oo), K.) such that: 

(22) ^ m I?([0,T],R). 

Then: 

^ )t>o « ([ 0 ,T],M+) , 

where {i^t)t>o deterministic element of C ^[0,T], defined for all f GVh 

and all t G [0, T] by: 

{i'tJ) = {i^o^nf)-ti[ Tt-sfifs) ds + X [ E[rt_s/(L>)] ds. 

Jo Jo 

Proof. Let (xt)t>o ^ limit point of < (p")t>o f • hand, assumption 

CD implies that (1 r ^ converges in distribution to 0. On another hand, 

L t”o n£N* 

for all (j) G Cl the mappings 


and 


: 


^'2 : 


I?([ 0 ,oo),] 

i^t)t>0 

V ([0, oo), I 

(Xt)t>0 


> V ([0, oo), K. 
■ {4>i^t))t>o 

C ([0, oo),M) 
{foYsds) 


t>0 


are continuous, as well as 'h := il '2 o 'Ll, hence in view of II22II . the continuous 
mapping theorem entails that: 


</'(^i(^")) ds 


f (l){rs)ds'\ inX>([ 0 ,r], 

Jo J t >0 


Consequently, for all t G [0,T], all cf) G Cl'. 

(23) {Xt,(l)) = {vlA) - i {Xs,(l)')ds- fj, f ())(fs) ds +AfE [(()(5)] . 

Jo Jo 

In particular, the latter is true for all (j) G S: (|2HJ is the integrated transport 
equation (E(pQ,g, 1)), where g is defined by: for all (j) G S, 

{gt,4>) ■=-g[ (j) (fs) ds + xm [(i){D)]. 

Jo 

According to Theorem^ the only solution of 12811 is given for all t G [0,T] and all 
4> G S by: 

ixt, (j)) = {vo,Tt 4>) + {gt^cf) - [ {gs,Tt-s (f') ds 

Jo 

= {pQ,Tt(j)) - g f 4> ifs) ds + AtE [4>iD)] 

Jo 

+ gf f {Tt-s4>'ifu))duds-x[ sE (/)'(£))] ds 
Jo Jo Jo 
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The limit point is therefore unique in V equal to since S is 

a separating class oi . □ 

7. Applications 


M/M/l/l+D-EDF case. We hereafter apply Theorem 0 to determine the fluid 
limit of the M/M/l/l+GI-EDF system in which the time credits of the customers 
are deterministic. We verify in particular that assumptions m and lp|l are sat¬ 
isfied in this case, and specify the form of the limit. 

We therefore consider a sequence of M/M/l/l+D-EDF systems, for which we 
make the following assumptions: 


Hypothesis 2 (Basic Assumptions for a M/M/l/l+D-EDF system). • For 
all n € N*, there are initially n +1 customers in the buffer, all of them with 
time credit nd, where d > 0 (that is, the n customers have their deadline 
at time nd), 

• for all n € N*, A" is the intensity of the Poisson process of arrivals, where 

A" —> A > 0, 

n^oc 

• for all n G N* the customers require service durations exponentially dis¬ 
tributed, of parameter pP, satisfying: 

/i" —> p , where {d)~^ < p < \, 

n—^oo 

• for all n € N*, the initial time credit of any customer is deterministic, given 
by d", where dP jn —> d. 


This queueing system is described by the profile process associated 

to {nSnd, , p^, d'^), which keeps track of all the residual time credits of all the 
customers waiting in the buffer are already lost. The notations are those of the 
preceding sections, with superscripts for “deterministic”. 


Lemma 1. For all x < p 


P [f-- < a:] 


0 . 


Proof. For the event 

f-n,D f n.D ^ 1 

{tq’ <a:| = |To’ <nx} 

to occur, the n + 1 customers initially present in the buffer must have all entered 
the service before time nx, since they couldn’t have been eliminated before nx < 
n/p < nd. Therefore the first n customers among them must have completed their 
service before nx, or in other words: 

P [fo’° < x] = P [t”’° < na:] < P [S'"^ > n ], 

where (5'r)t>o denotes a Poisson process of intensity /i” (the server works without 
interruption at least until he has completed the services of these customers). Hence, 
has the same distribution as the sum of n independent r.v (P",f = l,...,n), 
Poisson distributed of parameter p^x. Hence, denoting {Pi,i = 1, ...,n), a family 
of n independent r.v. Poisson distributed of parameter p, 


P < a:] < P 






1 

71 


2=1 




0 
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according to the weak law of large numbers. 


Lemma 2. For all C > 0, 


P r - n, D ^ D 

K <^^0 ' 


-*,D pd-fi ^ 


Proof. We have: 


(24) P - ^] = P [{Wo’° < - <} n {nWg’^ - < < Tg ’°}] 

+ P [{wo’° < Tg”’°} n {nWo’° - < > Tg”’°}] 

I Tj r f n.D ^ n,D ^ -*,D <^1 f -*,D > . n,D11 

+ P|_|rQ’ <U)o < nCJQ ~ > '^o jJ 

^ T> r ^ / -*,D . n,D1 I T> r n,D ^ n.D ^ -=k,D 

< F < (ntjQ — n^) A Tq J + F [Tq S ~ • 


Let us denote for alH > 0: 


On the event 


A/"" := number of customers arrived up to time t, 
having deadline before t 
5" := number of services completed up to t. 


( n.D ^ / -*,D . n,D'\ 

{Wg’ < (nWg’ -OATg’ }, 


there is at least one loss and no idle time before time (na)g’° — n^) A Tq’^ . For 
this to occur, since the service discipline amounts to FIFO, there must be at the 
first time of loss, say nt > nd, the number of customers initially in the system or 
arrived up to nt, of a priority higher to the priority of the customer who is lost 
at nt (there are n + 1 + such customers) must be greater than the number of 
services initiated up to nt (i.e., + 1). We can therefore write that: 


(25) P Wg’° < (nwg’^ — n^) A Tq’^ 


< P sup (Afnt + n + 1) — (5”( + 1) > 0 

d<t<(<:ug’^-0/\fa'^ 


= P sup AC - > -n 

Since there has been no idle time in [O, n(ujQ’^ ~ ?) A C’°], for all t in this interval, 
has the same distribution as S'”^, where (S'”)(>g is a Poisson process of intensity 
C- On the other hand, the process of arrivals (Ai")j>g marked by the initial time 
credits of the customers {L?"}ieN* being a two-dimensional Poisson process, it is 
easily checked that the process (A/)" — A”)t>o is an C”"ninrtingale, where 

A” = A” (t - nd)+ . 

Thus the process defined for all t by 

M7 ■.= N:,-Sl,-{\l,-p-nt) 
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is a -martingale. Hence, with (EH: 


(26) P [wq< (nwo’^ - O 


< P 


sup TW" > inf — A"^) — ^ 


sup 7Wr>n{A’^d-l-(A"-/r")(u;;’°-e)} 
4 




= P 


< 




rE 


^0 -?J 


!{A"d-l-(A"-^") (wo’^-e)} 

■ ,,.2 {»^" K-° -(-d)+ /.”n (oJ-> - ()} 

(n(A - 


using successively Doob’s inequality and the fact that (A — m)wo’° = Ad — 1. Now, 
clearly 

(27) P l_7-g ’ < oJq < nujQ — n^\ 


< P 


2/i 


on the event 


there exists t G 


2 fj, 

1 -*,D ^ 


^ n,D ^ n.D ^ —*,D /- 

< Tq S ^0 S ^^0 ~ 

^ ^ n,D ^ n,D ^ -*,D ^ i 

777 < T-g ’ < Wg ’ < nWg’ - , 


n,D ^ 
■^0 < 


2/x 


2 j 2 ,uuQ — (the first one) such that the buffer is empty at nt, 
but there has been no loss before nt. For this event to occur, there must be up to 
nt, the same number of customers entered (i.e., fV"^) as of services initiated (i.e., 
-I- 1). Thus, remarking that the process defined for all t by: 

:= iV- - 5” - (A”nt - ^"nf) 

is a C/"-martingale, and that for all t < equals in distribution. 


n 

2/7 


^ n,D ^ n.D ^ -*,D ^ 

— < Tq <^o — ~ 


< P 


< P 


sup 


sup — M" > —1 -f n(A" — — 

D t 2/7 


< 


4 


(„(An-/7»)^-l)' 

16 


rE 


< M” >_.,D , 

‘^0 -< 


(n(p - 1) - 2)' 


{(A" + /7")r7(:v*’°-0} ^ 0 , 
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which, together with Lemmas applied to a; = 1/(2^) and i?7ll yields: 

(28) P < Wq — n^l — s- 0. 

We conclude by substituting and in (E3. 

Proposition 2. Assume Hypothesis 1 and 2 holds, then 

where is the deterministic element o/C ([0, oo), K) defined for all t by: 

(29) . 


□ 


Proof. Fix ^ > 0, assuming without loss of generality that: 

(30) 

which implies that: 

(31) /i-i + pC < - C- 


d — p, ^ 1 

C< - 

P P 


If p < 3, assume in addition to pijl that: 


(32) 




1 


(4- p)p' 

Let us first focus on the interval of time [0,Wo’° — ^]. We have: 


(33) P 


sup \ti{v^' ) - ^ 




< P 


sup |ti(p°)-nrj^l > < 

0<t<| (lIio’°-J)AOo 


, T> r-=t^.D -n,D \ /-I I T> r ^ -*,D 

+ P [wq- - Wq ’ > 4 j + P [tq ’ < nwp’ - . 


Let us denote for all y G M and t > 0: 

^t,y •= number of customers arrived up to time t, having deadline before t + y. 
On the one hand, on the event 


I sup {ti{v^t) - nrf) > < > , 

|^0<i<|(ajQ’ —^)AajQ’ ATq ’ | J 

there exists t € [O, {(a)g’° ~ ?) A u>q’° A such that ti(i/lf^) — nfjr’ > n^. Hence, 

since there has been no loss until nt, the number of services initiated up to nt (i.e., 
iS”( + 1, which equals +1 in distribution) is larger than the number of customers 
initially present, or arrived up to time nt, having deadline before nt + nrf + 
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On the other hand, on the event 

(34) < sup {nf° - > <> , 

|^0<i<|(ajQ’ —^)AajQ’ ATq ’ | J 

there exists t S [O, (a)g’° ~ C) ^ such that < nfjr’ — n^. But for 

all t < (^“^) + p^, every customer initially in the system, or arrived before nt has a 
deadline at, or posterior to, nd, and hence a residual time credit at nt larger or equal 
to nd — nt > nf^ — n^. Therefore, for the event El to occur, there must exist an 
instant t S + p^, (u)q’° — 0 ^ ’°] ®uch that ti{n^^) < nf^ — n^. Since 

there has been no idle time on the interval [0,a;g’° A and since the discipline 

amounts to FIFO, this implies that + 1 (equal in distribution to + 1) is 
less than the number of customers arrived up to time nt, having deadline before 
nt + nrf — (that is, + n + 1). Consequently: 


(35) P 


sup \ti{v^t) - nrt \ > 

0<t<{ 


< p 


sup {tiiv°^)-nrf)>ni 

0 <t<| (iIIq’^-OAOo’^ATo’°} 


sup 

0<t<{K''^-OAcS"'UAr;'U} 


(nr° - ti{v°t)) > < 


< p 


sup (5;)* - V” > n 


0 <t<ag’^-i 


+ p 


sup 


(Ar..„(,p-fl - s;.) £ -» 


It is easily checked, that for all n and all z G R, 

is a t/”-martingale, where for alH > 0 : 


D+.) - ^m) 


t >0 


Kz = {t- {d- (rf + z)) 


r- 


Thus, the processes 


M-Pt) ■■= SZ - AJ,.„(,p+j) - (p”»( - AJj) 


»?(*) '= A:i,„(,p-£) - s?., - {k-( - 


and 
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are Qp— martingales. Hence, (EU becomes: 


(36) P 


sup - ur°| > < 


< P 


sup MP{t) > n + inf (a"£ — 

--•■o - ^ 


0 <t<Qp^-^ 

+ P 


sup Mp{t) > —n + 




inf (- XP A 


Ki \ „ I 2^7- _ 2u + A"f - 


On the one hand, since: 

n {A”(t A ^) —/r”t} iftG[0,/r“^] 

_,_1 _i_ _ ,,’^+1 

A A 

it follows from Hypothesis [2l that for a sufficiently large n, 

^A"£ — > n + n — I )■ >0. 


^ + A"C-/r"f} if f e [/r -^] , 


n+ ml 
o<t<ap^-i 


On the other hand, since for all t G ^ + p^, - ?] > 

- ~Xp_^ = nifi^t- -^t + — + A"^ 


for a sufficiently large n. 


-n + 


inf 


(p"nt - 


^'^nt — A”_£ ) > —n + n-(^ + l)’> 0 . 


p~^+pi<t<^o' “5 
Thus from for a sufficiently large n. 


(37) P 


< 


sup |ti(p°t) - nr°| > < 

0<i<{(io’°-€)A<I>J’°ATo”’D} 

16 




E 


{.M^ ((^*-° _ e) }^1 + E [{.M^ } 


16 


I A"., 


,^. + A” D +2p"n(w*’°-e)| -^0 
(nA^) 1 I'^o’ “f)A (“o’ } n^oo 

using successively Tchebitchef and Doob’s inequalities. 

Consider now the term: 

(38) P — n^] < P [wq< nwg’^ — n^] 

+ P [{tq< nuxl'^ — n > nwg’'^ — n^}] 


f n,D ^ -*,D r —*,D >1 

Itq <nujQ > ncjQ 


On the event 
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there exists an instant, say t G [0, such that the buffer is empty at nt, but 

there has been no loss before nt. Applying the same arguments that led to 
yields: 

P n {wo’° > - <}] ^ 0, 

^ ^ n—>'CXD 

which together with Lemma El and EHl) yields: 

(39) P \tq'^ < nu}l'^ — —> 0. 

^ n—^-<50 

Finally, using m together with and Lemma El in lp| : 


(40) 


sup 


I , / -n,D \ -D I 




Let us now consider the interval of time [ia>q’° — oo). We have: 


(41) 


P 


sup |ti(u”’°) - r°| > 2^ 


< P 


sup 






— nrf) > 2n^ 




/ \ 


+ p 

inf fi(u°)<n 

sup - 2^ 



t>i^o -i 




First, it is easily seen that 


sup 


t>U 




p-1 




Therefore: 





/ \ 


(42) 

P 

inf hi^^nt) < n 

sup ff-2^ 








On another hand, define the following event: 


:= < for all t > Wq’ — some customers arrive before n{t — ^), 


with deadline in [n{t + + ^), n(t + + 2^)] 
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We have: 


(43) P 


sup (ii(p°)-r°) > 2< 




< P 


sup (ii(p°)-nri°) > 2< 




n sup |ti(p°t) - nrf\ < nO n 




sup |ti(p”’°)-r°| > ^ 


0 < i < cS *’°-4 


P [{s^Y]. 


On the event 


sup (ti(p°) - nr°) > 2< 




n ^ sup \tiiv°t) - nrY < < ^ n EY 


there exists an instant t > a)g’° — ^ (the first one), such that there is no customer 
in the system at nt having deadline between n{t + rf + C) and n{t -\- rf + 2^). 
Consider the customer who would have had the smallest time credit at nt if we 
was still present in the system at this instant, among those arrived before n{t — ^), 
with deadline in [n{t + ff + ^), n{t rf + 2^)] . For all s < t — ^ such that this 
customer has already entered the system at ns, denote by Rg, the time credit of 
this customer at ns. We have: 


Rs G 



s) + nr° + n^, n(t 


s) + nr° + 2n^ 


In particular, it is easily seen with the form of that: 

Rs > n{t - s) + nf° + n^ > nf° + 

(the last inequality is true since s < Wq’° — ^ and 

sup |<i(n”;'^)-nf°| < nO- 

Thus, none of the customers arrived before n{t — ^), with deadline in 

[n{t + fP + ^), n{t + fP + 2^)] have been served before n{t — ^), since none of them 
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has ever been prioritary. Hence: 


(44) P 


sup (ii(u°)-nr°) > 2< 


L I 


n 


sup |ti(u°t) - nr°\ < < > n 




< P 


< P 


All of the customers arrived before n{t — ^), 

with deadline in [n{t + rf + ^), n{t -\-rf + 2^)] 

have entered service between n{t — ^) and 


sup I - f 




K(t-an(rP_^ + ^() ' ^ - 0 


where denotes the value at of a Poisson process of intensity The following 
is a tj"— martingale: 


:= - <! AT' 




-Af” 




|^”n^ - 2£+i^ - £+1, 


— 


)}■ 


and for all € L‘^. Then, it is easily seen that the last term of Ij44|l is equal 

to: 


. ' S 




But for all t > Wq’^ — 


\n \n 




2p + 1 






= A"C, 
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using m and Consequently, for a sufficiently large n: 


(45) P 


sup -nrt°) > 2< 




n 


sup nrf I < nO n 




< P 


sup A4"j > n (A" —/r") ^ 




< 


(n(A” -M”)0 

using again Doob inequality. Remark, that we also proved that: 
(46) P[(£-^")“] 

= P 


{A"n^ + Ai”n^} —> 0, 


< P 




sup +2^^) 


■^n(t-{).n(rP_^ + £±i{) j f” “ 


0 . 


Therefore, using (jm, (gni and M in o yields to: 


sup (ti(r'°)-r°) >2< 


Together with m in m, this entails: 


sup |ti(i^t”’ ) -'r°| > 2C 


This implies in turns, together with that: 


p 

sup ti(i/"’°) - > 3^ 

< p 

sup ti(!/"’°)-r° 


_t>0 




+ P 


sup |ti(i/"’'^) - r°| > 2^ 


—>0. 


We conclude using |Rohflflj . 


□ 
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We can therefore conclude with the next result, which yields then convergence 
in distribution of a normalized sequence of profile processes of a M/M/l/l+D-EDF 
queue to an explicit fluid limit: 


Theorem 4 (Fluid limit of the M/M/l/l+D-EDF queue). 




) 


t>0 


in V 




where for allt>0 and all f GVi,: 

ds + X f {d — s) ds, 

Jo 

i^F)t>Q being defined by 


F‘^,f)=fid-t)-nf /(rf + s-t) 
Jo 


Proof. Let us verify the assumptions of Theorem |31 for any given T > 0. First, it is 
straightforward, that Hypothesis Hare satisfied in view of Hypothesis |2l for D = d, 
a.s., and for ^ the Dirac measure at d. 

Then, remark, that in (HKll w 6 prov 6 cl cls cl nicittBr of fcict thcit for till ^ 0 

satisfying PI and (ESI, 


There exists t > uJq’^ — such that all of the customers arrived 
before n{t — ^), with deadline in [n{t + ff + n{t + ff + 2 ^)] 

have been served before nt 




0 . 


Denoting this previous event it follows that: 

P [ro< T] < P - C] + P [Al + P 

using pji and pjl as well: EU is satisfied. Finally, (| 22 j is verified in view of 
Proposition [21 for (d°)i>o defined in 12311 . We therefore can apply Theorem (31 for 
all r > 0 , which completes the proof. □ 

We can in particular approximate the queue length and loss processes, by apply¬ 
ing the profile process of the queue to simple rcll functions. First, in view of m, 

the normalized queue length process can be asymptotically approximated 

by the process deflned for all t > 0 by: 

, 1r* ) = -M / Ik* {fF + s-t) ds + X / 1r. {d - s) ds 
Jo Jo 

{1 “t” ^ f pd-ij.-^ j “t“ Xdl. f j ■ 

Similarly, in view of dJ the normalized loss process (P")j>q can be approximated 
by the process deflned for alH > 0 by: 

1r_) = -d f 1r- {rf + s - t) ds + X [ 1r_ {d - s) ds 
Jo Jo 

= (1 + X{t -d)- fit) • 


/ -*,D 

w 
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M/GI/ oo system. Remark, that we can also apply Theorem|2|to obtain the fluid 
limit of a pure delay M/GI/oo system. For all n S N*, consider a pure delay (PD) 
M/GI/oo system: customers arrive according to a Poisson process of intensity A", 
requesting service durations, i.i.d. of the distribution of a", to an inflnite reservoir 
of servers, assuming that: 


A” 


A, E 


a 

n 


E q; ,-> a. 

n 


Each customer is hence immediately attended upon arrival, and remains in the 
system for the duration of his service. Such a system can be described by a profile 
process, keeping now track of the remaining processing times of the customers in 
service (positive atoms) and the elapsed times since departure of the already served 
customers (negative atoms): this is the profile process associated to 

(pqA", 0, q;") , where Pq is the profile of the service durations of the customers 
initially in the system (which we assume to be an n-sample of the distribution of 
a”). The service durations replace the time credits of the queue with impatient 
customers, and consequently the analogous service rate becomes null in this case. 
We normalize this process as above, writing for all t > 0 and all Borel set *B: 


-n,PD 




m = 


n,PD 


(nQ5) 


The fact that /i^ is zero allows us to skip conditions eu and l(23) in Theorem 0 
whose application becomes straightforward: 

{i^t ^ ([0,cx)),At+) , 

where deterministic element of defined for all t > 0 and 

all / e Df, by: ■ 

ct 


-*,PD 




> /) = E [nfia)] + A / E [Tsf{a)] ds. 
Jo 


Hence, as above, we can asymptotically approximate the normalized congestion 


process (number of customers in service), given for all n e N* by 
by the process defined for all t > 0 by: 


-n,PD 


t>0 


A 1 

, J-IE 


) = P [a > t] + A / P [a > s] ds. 
Jo 


Remark, that in the special case where a" is exponentially distributed of parameter 
/x" ^ this process becomes for alH > 0 : 

e-^‘+ p (1 - , 

which is the fluid limit obtained by Borovkov in |Bor M- We can also approximate 

by the 


the normalized workload process, given for all n by ^( 
process given by: 


-n,PD T-i 


t>0 


/-=*'.PD r-i 


) = E [(a - t)+] + X [ E [(a - s)+] ds. 
Jo 
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Finally, the normalized process counting the already served customers can be ap¬ 
proximated by the process defined for alH > 0 by: 

1r_) = P [a < f] -|- a / P [a < s] ds. 

Jo 
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